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Abstract 

A general framework for deconvoluting the effects of energy averaging on 
chargcd-particlc reaction measurements is presented. There are many poten- 
tially correct approaches to the problem; the relative merits of some of are 
discussed. These deconvolution methods are applied to recent ^^C(a,7)^®0 
measurements . 
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1. Introduction 

Cross section measurements in nuclear physics in general involve an average 
over a range of energies. The importance of these convolutions depends on how 
strongly the cross section varies with energy as well as experimental conditions. 
In the case of experiments with charged-particle beams, which is the intended 
focus of this work, the most important experimental effect is usually beam 
energy loss in the target. If these effects are significant, the data must be 
corrected for them before they can be used (e.g., for astrophysical calculations, 
or comparison to other experiments or theoretical calculations). While these 
effects have been known by workers in the field for many decades, a general 
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treatment has not been given and some confusion remains in the hterature. In 
a recent paper, Lemut [1^ has advocated a particular approach for correcting 
experimental data for energy averaging. It is, however, but one of many valid 
approaches, and, in some circumstances it has some distinct disadvantages. In 
this paper, we present a general framework for deconvoluting the effects of 
energy averaging. We emphasize that there can be many correct approaches 
and discuss some of their advantages and disadvantages. 

The ^^C(a,7)^^0 reaction is a very important process in nuclear astro- 
physics. Differential cross section data for this reaction have been recently 
published by Assungao et al. ^2j. This data set provides an opportunity to 
demonstrate energy deconvolution, which also renders the data in a form that 
is more useful for further future analysis. 

2. Energy Deconvolution 

2.1. Statement of the Problem 

For an ideal experiment, with a mono-isotopic target of uniform areal density 
and with perfect energy resolution, the reaction yield Y (number of reactions 
per incident particle) is related to the cross section a via 



where n is the number density of target atoms. Ax is the linear target thickness, 
and E is the beam energjj^ In a real experiment, a variety of additional physical 
effects must often be taken into consideration. For experiments performed with 
charged-particle beams, which are the focus of this work, beam energy loss in 
the target is an important effect. Considering energy loss in the target, the 
relationship between yield and cross section becomes 



Y[E) = (j[E)nI^x, 



(1) 



Y{Eq) = / (j{E{x))ndx 



L 



dE, 



(2) 



^We utilize laboratory energies throughout this paper, unless otherwise indicated. 
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where Eq is the incident beam energy, x measures the hnear depth in the target, 
AE is the energy loss in the target, and e = — is the stopping power. The 
stopping power is energy-dependent, and AE must generally be determined by 
a numerical procedure. For example, the variation of the energy with depth in 
the target can be determined by numerically solving the differential equation 

^(x)^-neiE), (3) 

subject to E{0) = Eq, and then AE — Eq—E{Ax). The important point that we 
are concerned with here is that now the yield is related to an energy convolution 
of the cross section. Energy loss in the target is usually the most important effect 
contributing to the overall energy resolution in an experiment, but additional 
effects such as energy straggling, the incident beam energy resolution, energy- 
dependent detector efficiency, and target composition may also be important. 
These additional effects lead to a more complicated form of Eq. ^ that in some 
cases is most practically implemented via a Monte Carlo simulation. Since the 
nature of these effects depends upon the details of the particular experiment, 
we will not consider them further here. It should be noted that considerations 
discussed below can be generalized in a straightforward manner to these more 
complicated energy convolutions. Equation ([2| can also be applied to differential 
cross section measurements, in which case Y becomes the yield per steradian 
and cr must be replaced by the differential cross section. 

The question we want to address here is: how does one convert Y{Eo), the 
yield measured in an experiment, to a cross section? In other words, how does 
one deconvolute the effects of the energy averaging process? Our primary inter- 
est is to compare the relative merits of various mathematical approaches that 
have been posed in the literature. In order to carry this out, it is necessary that 
the energy dependence of the cross section and the physical processes generat- 
ing the energy convolution be well understood. We will, for now, assume these 
assumptions are fulfilled. The important issue of error analysis arising from 
uncertainties in the experimental data and the energy convolution is dependent 
upon the details of the particular experiment and will be discussed briefiy at 
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the end of this section. The effects of angular convolution, resulting from the 
acceptance of the detection system, are not considered in the analysis presented 
in this section. 

These deconvolution corrections are expected to be most important when the 
cross section is highly energy dependent (e.g., charged-particle reactions below 
the Coulomb barrier or resonances) and/or when when thicker targets are uti- 
lized. It should also be noted that the methods discussed in this paper may not 
be the most appropriate for analyzing data from narrow resonances (with widths 
significantly smaller than AE); in those cases one generally determines the reso- 
nance parameters (e.g., resonance strengths and resonance energies, rather than 
cross sections) directly from the measured yields. 

2.2. The Effective Energy Approach 

One approach is to to utilize the ideal formula, Eq. Q, to define the exper- 
imental cross section 

YjEp) 
n Ax 



<Tcxp(i?cff) = (4) 



where Ecs is the effective energy^ defined in this work to be given implicitly by 

'j^^E. (5) 



In Eq. ([5]), the assumed theoretical function for (7{E) is used on both sides 
of the expression; the solution for E^s thus requires inverting the cross section 
function. Since this approach ensures that the measured yield and experimental 
cross section are consistent with Eq. ([2]), this procedure can be considered a 
faithful (mathematically-correct) deconvolution. As far as we are aware, this 
approach was first described by Dyer and Barnes in 1974 [3], and has been 
used by several authors subsequently (see, e.g., Refs. [11 15]). It has also been 
described and advocated for in a recent paper by Lemut [T]. 
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2.3. A General Class oj Approaches 

A more general class of deconvolution procedures can be derived by consid- 
ering adjustments to the experimental cross sections in addition to the energy: 

Here, E is the energy assigned to the experimental measurement and the cor- 
rection factor / is given by 

^EUE^{E)[<E)]-'dE 

The factor / depends upon the assumed energy dependence of <j{E) and again 
the procedure is consistent with Eq. ([2]) , ensuring a faithful deconvolution. The 
method for determining E is not specified at this stage, but one expects rea- 
sonable procedures to provide E within the range of energies in the target, i.e., 
Eq — Ai? < E < Eq. One can consider the effective energy approach described 



in Subsec. 2.2 to be a special case where E is adjusted to make / = 1. 

An intuitive choice for E is to define it to be the mean or cross-section 
weighted energy: 

^_ lE:-AEEa{E) [ejEr^dE 
lE:-AE'^(E)[eiE)]-^dE ■ 
To the best of our knowledge, this procedure was first used by Wrean et al. 
[B] . A similar approach is described by Dwarakanath 7 . This choice ensures in 
some sense that the energy assigned to an experimental cross section measure- 
ment reflects the most important reaction energies in the target. An additional 
utility of this method is that the mean energy can be measured directly in some 
cases, via the detection of capture gamma rays [5HlO|. This approach was mis- 
characterized in Rcf. [1 , as the correction factor / was left out of that author's 
description. 

Another possible choice for E is to define it be the median energy £"1^, which 
is defined implicitly: 

''^'^^dE^2r^-j$dE. (9) 

Eo-AE 



e{E) e{E) 
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This definition has been given in Refs. [TTHT5] . but these authors have wrongly 
imphed that this median energy can be used together with Eq. Q to faithfuUy 
deconvolute yields into cross sections. From the discussion above, it is clear that 
the median energy is of no particular significance in the deconvolution process 
and that Eq. ([6]), which includes the deconvolution factor /, must be used to 
faithfully deconvolute yields into cross sections. This approach was recently 
used by Makii ct al. 14 for the analysis of ^^C(a,7)^®0 measurements. In 
this case, the deconvolution factor is significant and is taken into account, with 
/,„ « 0.77 for both energy-target combinations reported. 

2.4- Analytic Approximations 

We will briefiy consider here some analytic approximations to the deconvolu- 
tion approaches given above. We ignore the energy dependence of the stopping 
power and assume the energy dependence of the cross section to be at most 
quadratic. Expanding the quadratic around the energy at the center of the 
target, = Eq — AE/2, we have 

where at = cr{E}i) and the first and second derivatives are also evaluated at Eh- 
It is also useful to define 

_ , \ and /?n = f - 

cJh \AE 

For the case of the effective energy, Eq. ^ leads to a quadratic equation: 



/dVN /d. 
\dE'^ ) ' \dE 



11) 



i?2e' + 2e- ^i?2(A£;)' = 0, (12) 
where E^fi = Eq ~ AE/2 + e. The resulting effective energy is 

E.,.^;.-> + £±MMW^^, (13) 

Z 1X2 

where we have adopted the solution which approaches Eq in the AE — limit. 
This point will be discussed further below. It is interesting to note that if 
i?2 = 0, i.e., if the cross section depends only linearly upon the energy, then 
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using EcS = Eo — AE/2 provides an exact deconvolution (the energy dependence 
of the stopping power is also ignored). 

For the case of the mean energy, we set i?2 = 0, due to the more comphcated 
integrations required. The results are 

E ^ Eo - ^AE + ^RiiAE)^ (14) 



and 



r' = l + ^RliAEf. (15) 



For the the case of the median energy, we also take i?2 — 0, resulting in 

^„ = go-lA^+ " + "''^f/^''''-' (16) 
Z Hi 



and 



J 111 



l + lRliAEf 



1/2 



(17) 



Equation ( 16 1 is mathematically equivalent to Eq. (5.55) of Ref. [T^], but, again. 



that reference did not consider the need for the deconvolution factor fm- 

2.5. An Example Case 

A good example of the application of these methods is provided by the 
^Be(a,n) reaction, which was measured for 0.23 < Ea < 2.7 MeV by Wrean 
et al. [6]. Over this energy range, the cross section includes the effects of the 
Coulomb barrier as well as narrow and broad resonances. The cross section 
varies by eight orders of magnitude over this range and is shown in Fig. [l] 
Three targets were utilized for the measurements reported in Ref. [S] , consisting 
of 5.08 X 10^^ (target 1), 1.23 x 10^^ (target 2), and 3.23 x 10^* (target 3) 
atoms/cm^ of ^Be, respectively. These data were analyzed using the mean 
energy method described above, with the energy dependence of the detector 
efficiency also taken into consideration. Consistent results for the cross section 
were deduced near the 0.618-MeV resonance using all three targets, and for 
energies down to Ea ~ 0.22 MeV using targets 2 and 3. 

Using the energy dependence of the cross section defined by the parameters 
in Table I of Ref. [6 and also the stopping power from that reference, we have 
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Figure 1: The ^Be(o,?i) total cross section, as a function of laboratory a energy, calculated 
from the parameters given in Table I of Ref. [6]. 

calculated the mean energies and corresponding correction factors for these tar- 
gets using Eqs. ([7]) and (|8|; we have ignored the beam-energy dependence of 
the detection that was also taken into consideration by Ref. [6j. The result- 
ing correction factors are shown in Fig. [2] for low energies. As expected, the 
correction deviates more from unity for thicker targets and/or when there is a 
strong energy dependence of the cross section due to resonances or the Coulomb 
barrier. The results using the median energy are qualitatively similar and will 
not be discussed further. For the effective energy approach, / = 1 and all the 
deconvolution is contained in the effective energy. 

2.6. Discussion 

It is found that there are many possible approaches to the energy deconvolu- 
tion that are mathematically exact. However, in practice, these different choices 
may lead to different results, due to the fact that the energy dependence of the 
cross section and the nature of the energy convolution are not exactly known. 
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Figure 2: The correction factor / for target 1 (solid curve), target 2 (dotted curve), and 
target 3 (dashed curve). 



It may also be the case in practice that approximate deconvolutions, such as 
ignoring the correction factor /, may provide perfectly acceptable results. 

In the preliminary stages of the analysis of the data presented in Ref. [6] , the 
effective energy approach described above was attempted. Two shortcomings 
of this method were encountered. The first problem is related to the fact that 
it requires inverting the cross section function. If the cross section has a local 
maxima or minima (e.g., at the peaks of or between resonances), this inversion 
is ambiguous. A criterion for selecting the solution must then be adopted; 
choosing the one nearest Eq — /S.E/2 is one reasonable convention. In any case, 
one finds that EcS_{Eo) is a discontinuous function and that certain ranges of 
E^fi cannot be reached for any Sq. Effectively, data points are pushed away in 
energy from cross section maxima or minima. The analytic approximation for 



the effective energy given by Eq. (12 1 is a quadratic equation and provides an 
example of this ambiguity. If the cross section has a local extrema at the center 
of the target, i.e., if {Aa/dE)E^ = 0, then ^eff = E^- l\E /2±/\E / ^/V2., and the 
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Figure 3: The dependence of the mean energy (soUd curve) and effective energy (dashed curve) 
on the target thickness is shown in the upper panel, for ®Be(a,n) and an incident a energy 
of 0.24 MeV. The lower panel shows the corresponding fractions of the yield resulting from 
reaction energies above the mean energy (solid curve) and above the effective energy (dashed 
curve) . 
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choice of roots in ambiguous. In the case of resonances (maxima) this behavior 
was not deemed desirable because the presence of the resonance was largely 
responsible for the measured yield and one expects reasonable procedures to 
assign an energy to a measurement that reflects the reaction energies that gave 
rise to the measured yield. 

Another shortcoming noted is that for thicker targets at energies well below 
the Coulomb barrier, is located at an energy significantly less than the 
energies where the majority of the yield is generated. In Fig. [3j the mean 
and effective energies are plotted versus the target thickness, for the ^Be(a, n) 
reaction with an indicant a energy of 0.24 MeV (approximately the lowest energy 
measured in Ref. [6 ). It is seen that the mean energy approaches a constant 
value once a thickness of « 2 x 10^*/cm^ has been reached, as additional target 
thickness does not contribute to the integrals in Eq. ([s]). However, the effective 
energy continues to decrease as the target thickness is increased, as additional 
target thickness does continue to contribute to the integral on the left-hand side 
of Eq. ([5]) . The lower panel of Fig. [s] shows the fraction of the yield resulting 
from reaction energies above mean and effective energies. It is seen that for 
the mean energy the fraction remains nearly one half, while for the case of the 
effective energy it approaches one. This behavior of the effective energy is not 
desirable, as the effective energy provides a misleading picture of what reaction 
energies were actually measured. In this situation the effective energy approach 
will also be more prone to systematic errors resulting from imprecise knowledge 
of the energy dependence of the cross section. 

It was for the above reasons that the effective energy was not used by Wrean 
et al. j^. Another situation where the effective energy approach is less useful 
is when multiple yields are analyzed from a single bombardment, such as for 
cross sections to different final states or differential cross sections for different 
angles. The energy dependences of these cross sections may be quite differ- 
ent, e.g., because of resonances involving particular partial waves, resulting in 
different effective energies being assigned to the different measurements. The 
resulting data are much more difficult to present and further analyze: consider. 
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for example, an angular distribution where each angle corresponded to a differ- 
ent energy! This complication can be avoided via the mean or median energy 
approaches, using the total cross section to calculate the mean or median energy, 
and then applying the appropriate correction factors to the individual partial 
or differential cross sections. This approach is used below for the analysis of 
^^C{a,'y)^^0 angular distribution data. 

The above difficulties with the effective energy method are related to the 
fact that all of the deconvolution is performed by adjusting the effective energy. 
These difficulties can be avoided by adopting alternative choices for the energy 
and determining the experimental cross section via Eqs. ^ and Q. Both the 
mean and median energies, defined by Eqs. ([s]) and ([9]), are reasonable choices 
that assign energies to the experimental cross section that accurately reflect 
the energies that were measured. We are not aware of any cases where there 
would be significant differences between these two approaches. We have a slight 
preference for the mean energy, as it is easier to calculate than the other energy 
definitions discussed here. 

In all of this discussion, it has been assumed that the convolution kernel and 
the energy dependence of the cross section are known. In practice these inputs 
may not be precisely determined, and error analysis must be considered. For 
the cross section extraction, it has been found that using self-consistent iterative 
procedures to determine the energy dependence of the cross section work quite 
well in several cases; see, e.g., Lemut '!], Wrean et al. [B]. Uncertainties arising 
from the cross section assumptions and the possibility missing features in the 
cross section due to insufficient energy resolution must be considered. Partic- 
ular attention should be paid to the lowest-energy cross section data point (or 
points) , where the assumed energy dependence of the cross section is obviously 
not constrained by any lower-energy data points. In addition, the energy regions 
above strong resonances should be scrutinized for the possibility that a tail in 
the energy resolution function causes measurements above the resonances to re- 
ceive resonant contributions. The latter effect is important in the case discussed 
in the following section. 



12 



3. Application to ^'^C{a,'y)^^0 data 

Extensive angular distribution measurements of the ^^C(a,7o)^^0 reaction 
were performed at Stuttgart a few years ago using detectors from the EU- 
ROGAM array [21 HI] . Measurements of this type are necessarily taken with 
targets which give rise to significant energy averaging effects. The energy decon- 
volution is complicated by the fact that the energy dependence of the differential 
cross section depends upon the angle, particularly near the Ec.m. — 2.68-MeV 
resonance. In Refs. jUIH], a complete energy deconvolution was not attempted. 
The primary goal of the present analysis is to more fully take into account 
the effects of energy averaging and to put the data into a form that is useful 
for further analyses, including a simultaneous _R-matrix analysis of the angular 
distributions. New information about the interference structure involving the 
Ec.m. = 2.68-MeV resonance [TB] is taken into account. 

3.1. Cross Section Model 

In order to consider the effects of energy convolution on the differential cross 
section, a model for its energy dependence must be adopted. The differential 
cross section for ^^C(a,7o)^^0 can then be written as [51 [T7] 



da- _ 1 



<JEiil - P2) + ^7^2(1 + - yi"4) 



+6cos0j'^(Pi-P3) 



(18) 



where aEi and aE2 are the El and E2 ground-state cross sections, </> is the 
relative phase, and the angular dependence is contained in the Legendre poly- 
nomials Pi = Pi{cos9). The term cx yfoEio^ results from E1-E2 interference, 
and gives rise to an asymmetry around 90° in the differential cross section. It is 
often convenient to to use the S'-factors Sel = o'elE exTp{2TTri), where 77 is the 



Coulomb parameter. Equation ( 18 1 can also be written as 

4 

W(f)]^ with W(f)] = 1 + ' 

dO 47r 



'^'^ = ^Wie), with Wie) = l + J2aiPi, (19) 



1=1 
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Figure 4: The R-matrix parametrizations for Sei, Se2, and cos </> as a function of CM energy. 
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where ctq = <7ei + <yE2, 



6 cos (j) l(TElCrE2 



To V 5 



(20) 



oi = -as 



—7aEi + 5crE2 



, and 



(21) 



7(7o 

12CTB2 
7 (To 



(22) 



If the above theoretical expressions for the differential cross section are to be 
compared to experimental differential cross sections measured with a finite sohd 
angle, Pi should be replaced with QiPi, where the Qi arc the attenuation fac- 
tors [IS]. We will utilize Wq{9) and W{9) to indicate angular dependences with 
and without the attenuation coefficients, respectively. 

For computing the convolutions, we assume that ctei, crE2, and cos0 are 
given by _R-matrix parametrizations, as shown in Fig. [4j where for convenience 
iS* factors are shown in place of cross sections. The El parametrization is essen- 
tially the same as that described in Ref. !19j, except that more recent phase shift 
data [3D] has been fitted. The E2 parametrization is taken from Ref. [16^ (the 
solid curve in the lower panel of Fig. 2) . The parametrization for cos (f> is taken 
from the relative phase of the above El and E2 parametrizations. It should 
be noted that cos (f> is primarily defined by the P- and D-wave scattering phase 
shifts [2T] and that our parametrization accurately reproduces the recently- 
published phase shifts of Tischhauser et al. [20'. Many of the details of these 
parametrizations are unimportant for the convolutions discussed below. The 
essential ingredients are the roughly constant S'-factors below Ec.m. — 2 MeV, 
and the behavior of <jeiI<^E2 and cos0 near the narrow Ec.m. = 2.68-MeV 
resonance. This behavior is largely determined by the known resonance param- 
eters, but the nature of the interference was until recently not understood |21j . 
However, the experiment and analysis of Sayre et al. [TOJ have now conclusively 
determined how this resonance interferes with other El and E2 amplitudes. 

3.2. The Experimental Data 

Refs. [2 [TJ] report measurements at 25 bombarding energies that are de- 
scribed in Tables 1 and II of Assungao et al. [2j particular energies will be 
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referred to by measurement numbers 1-25 which are ordered by energy as given 
in these tables. Gamma-ray angular distributions were measured at nine angles 
and are given in Figs. D.1-9 of Fey |15 . These data are the efficiency-corrected 
numbers of gamma rays detected, scaled by 10^^; these quantities are defined 
here to be N^{9). The data have also been corrected for finite solid angle. When 
normalized by the factor 0.624(47r)QiV, where Q and N are given in Table II of 
Ref. [2], these data become differential cross sections in nb/sr. We have verified 



that these differential cross sections, when fitted with Eq. (18 1, reproduce the 
cross sections and their uncertainties given in Table I of Ref. 2 , for both the 
two-parameter and three-parameter approaches described therein. The phase 
(j) and its error was also well reproduced in the case of the three-parameter fit. 
The fit was determined by minimization, with the parameter errors indi- 
cating the range of parameter values consistent with an increase in by an 
amount < 1 above the minimum when the remaining free parameters are var- 
ied. The uncertainty in N is propagated after fitting the angular distributions. 
The uncertainty in center-of-mass (CM) energy, given in column 2 of Table 1 of 
Ref. |2], is propagated to the S factors when calculating them from the cross 
section. 

These fits also describe the reported angular distribution data via N~^{d) = 
NjoW{9). We have removed the correction for finite solid from the data by the 
transformation N^{9) -> N^{e) + N^o[Wq{9) - W{e)], where A^^o and ctei/ctei 
are taken from the two-parameter fit. This approach by construction ensures 
that the identical results will be produced when the modified N^{0) is fitted 
with the attenuation coefficients taken into account. This step was taken be- 
cause, in the simultaneous fits of multiple angular distributions, we feel it is 
not desirable to have already made corrections that are based on single-energy 
fits. It also makes it possible to investigate the effects of uncertainties in the Qi. 
These factors are given in Ref. [5], having been calculated from a simple for- 
mula, their Eq. (4.2). We have calculated the Qi independently, using a GEANt4 
simulation ^^l, taking into account the tapered geometry of the EUROGAM 
detectors [23l|24] and the angle-dependent detector distances (see Table 3.2 of 
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Ref. [H]). The changes in the extracted cross sections from using our Qi were 
generally found to be negligible. The one exception was for the four El cross 
sections underlying the narrow E2 resonance (measurements 18-21). However, 
these four El cross sections also have very large statistical uncertainties, ren- 
dering the changes to be of less importance. We have thus used the Qi factors 
given in Ref. [2 in all subsequent analysis. 

Although the "effective energy" is extensively discussed by Assungao et al. 
P], Refs. [HIH] do not explicitly state how it is defined. They presumably used 
what is called in this work the median energy, as that was utilized in a slightly 
earlier analysis by the same group; see Eqs. (3.19) and (3.20) of Ref. ps] . 

3.3. Center of Mass Motion 

The analysis given in Rcfs. |H[TS] did not consider the effects of CM motion. 
We have included it in our analysis, using Eq. (B9) given by Brune [51]. This 
approach leads to a modified form for the angular distribution: 

4 

WQp{e) = Pqp^ + J2aiPQpi{cos9), (23) 
1=1 

where 

Pq^o = l-H2/3QiPi, (24) 
PqPi = QiPi + ^^[{l + l)il + 2)Qi+iPi+,-l{l-l)Qi^,Pi.,], 

and f3 is the speed of the recoiling ^^O nucleus relative to the speed of light. We 
have utilized a value of 0.849 for the required attenuation coefficient. The 
inclusion of these effects makes some non-negligible changes, particularly for the 
extracted E2 cross sections with Ec.m. between 2.0 and 2.6 MeV, which are re- 
duced by 10-15%. The reduction is very similar to the case reported in Ref. pT| . 
The effects of CM motion have been included for all analysis described below, 
with f3 calculated for the median energies in the target reported by Ref. [2] . We 
have verified that this approximation (ignoring the energy dependence of /3 over 
the beam-energy loss in the target) is excellent. It should be noted that the 
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resulting formula for the differential cross section remains a linear combination 
of three terms proportional to asi, <^E2, and cos (t>^(JEi'^E2- With the assump- 
tion of constant /3, all of the energy dependence of the differential cross section 
is contained in these three quantities. 

3.4^. Energy Convolution Model 

The targets consisted of ^^C implanted into gold. The depth distribution can 

be described by g{y), where g is the ratio of ^-^C to (^^C + Au) and y measures 

the depth in units of (^^C + Au) atoms/cm^. Equation [2] becomes 

rva 

Y{E^)= / a{E(y))g{y)dy, (25) 



"'0 

where E^ is the incident a energy, yo is the maximum depth considered (where 
g{y) — >■ 0), and the variation of energy with depth is determined by solving the 
differential equation 
AE 

-^{y) = -5(y) ec{E) - [1 - g{y)] SAuiE), (26) 

subject to E{0) = Ea, where ec and eau are the stopping powers for helium 
ions in carbon and gold, respectively |26) . The total number of carbon atoms 
per unit area in the distribution g{y) is given by 

Ng^ / g{y)dy. (27) 
Jo 

We have neglected energy straggling, as investigations into this process indicated 
that this effect is small compared to changes associated with uncertainties in 
g{y). The calculated target-averaged differential cross section can thus be de- 
fined to be 

f";;^(i^)5(.)d.; (28) 



' dcr \ 1 ry° da 

recall that the corresponding experimental value was noted in Subsec. [3.2|to be 



da\ _ N^ie) 



(29) 

dr!/,,p 0.624(47r)giV 
We have thus allowed for the possibility that g{y) has a normalization that does 
not correspond to the experimentally-determined areal density of carbon, i.e., 
that N ^ Ng. 
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3. 5. Three-Parameter Fits Near the Narrow E2 Resonance 
I.Of 




18 2 

y (10 atoms/cm ) 



Figure 5: Depth profiles g{y) for the targets used near the narrow E2 resonance. The before 
(sohd curve) and after (dotted curve) profiles are taken from Fig. 8 of Ref. [2]. The dashed 
curve is discussed in the text. 



We begin by considering measurements 13-25 which are near the narrow 
E2 resonance and were conducted with very similar targets. For bombarding 
energies at and above the resonance energy, the effects of energy convolution 
are substantial. At this stage, we will proceed by comparing measurements 
of the target-averaged cross sections to calculated convolutions, as opposed to 
proceeding straight to deconvolutions. 

A three-parameter fit to the experimental target-averaged differential cross 
sections may now be performed using 



da 
dO 



1 

An 



00 



QI32) 



130 



7 



QI32 



-1-6 (cos ( 



[cEl){<yE2) 



1 - Pqp^) 



to yield experimental values for (aEi), {o'e2), and (cos( 



(30) 

The theoretical ex- 
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Figure 6: Results for the experimental and theoretical {a el) a-^d (cos</>) versus bombarding 
energy. The results from measurements 13-25 are shown as the circles and the theoretical 
calculations are shown as the curves, which correspond to the depth profiles shown in Fig. [5] 
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pectations for these quantities are: 

1 fy 

{<^el) = / '^EL9iy)dy (31) 

JVg Jo 

1 fy 

(C0S(/)) = ~ ,, / COS(?!'Vcr£;iCrE2g(y)dy, (32) 

NgyJ {oei){(JE2) Jq 

where L = 1, 2 and the energy dependences of the cross sections and cos (j) under 
the integrals have been suppressed. 

These measurements were carried out with targets containing approximately 
1.3 X 10^* ^^C/cm^. Example before-and-after depth profiles of such a target 
are given in Fig. 8 of Ref. and are also shown in Fig. [sj The results for the 
experimental and theoretical {(J el) and {cos(f>) are shown in Fig.joj In addition 
to the two depth profiles from Ref. ^ , we have considered a third profile that was 
adjusted to optimize the agreement with the experimental (cr_B2)- We find that 
{(7ei) is insensitive to the details of the depth profile for all energies considered. 
Below the resonance energy, {<te2) and {coscf)) are also insensitive to the details 
of the depth profile. However, at and above the resonance energy, (cos 0) and 
particularly {aE2) are very sensitive to the depth profile. 

The third depth profile, shown by the dashed curves in Figs. [5] and [6j includes 
a small tail which extends beyond a depth of 2.5 x 10^* atoms/cm^. It is this tail 
which allows the calculation to reproduce the two highest-energy data points. 
Note that the two depth profiles from Ref. P] do not allow the resonance to 
contribute to the {<te2) calculation at the two highest energies. We believe 
that is extremely likely that E2 contribution for these data points is in fact 
coming from the resonance and the tail of the depth profile. The alternative 
explanation would be that our model cross section is too low at these energy 
by a factor of ^ 5, which is very unlikely (see Ref. ll6J). In addition, the 
Rutherford backscattering technique, which was used to determine the depth 
profiles in Ref. [5], is not sensitive to a small tail. Further evidence for this 
interpretation is provided by the gamma-ray spectra given in Figs. C.24 and C.25 
of Ref. [15] that correspond to these two data points. The full-energy peaks from 
the '^^C(a, 7o)"'^^0 reaction at angles away from 90° are double-peaked. The 
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lower-energy portions of the peaks match the peak energies for the on-resonance 
measurements, an observation consistent with significant reaction yield coming 
from the resonance. 

We conclude that meaningful E2 cross section cannot be extracted for en- 
ergies at and above the resonance energy {Ea — 3.580 MeV and higher) due 
to the strong sensitivity to the depth profile. Similar conclusions were reached 
by Assungao et al. 0, except that they extracted cross sections and S factors 
for the two highest-energy measurements. For the reasons given in the previous 
paragraph, we do not believe these two points can be reliably analyzed. The 
El cross sections can separated from the uncertain E2 contribution using the 
above three-parameter fit procedure; the {(Jei) will be analyzed further below. 

3.6. Deconvolution of the Angular Distributions 




Figure 7: The depth profile g(y) for measurements 1-10. 



The depth profiles utilized in the analysis were defined as follows. For mea- 
surements 1-10, which utilized relatively thick targets, we have assumed the 
depth profile shown in Fig. uj In Ref. [2], median energies were determined 
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using the assumption of a constant S factor. We have determined the over all 
scale of these depth profiles by adjusting the maximum depth yo to reproduce 
the the median energies given by Ref. [2 , where we have, for this part of the 
calculation only, assumed a constant S factor in order to reproduce the methods 
of Ref. [2 ■ It should be noted that this procedure generates depth profiles with 
""^^C areal densities that agree with those given in Table II of Ref. |2j to better 
than 10%. For measurements 11-25, which utilized relatively thinner targets, 
we have assumed the depth profile shown by the dashed curve in Fig. [5| In 
these cases, the depth scale of each profile has been scaled by A^/1.30 so that it 
reproduces the experimental carbon areal density. 
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Figure 8: The calculated correction factors for 48° (circles) 90° (squares), and 130° (diamonds) 
versus measurement number, for measurements 1-17. 



In order to minimize the number of changes in the analysis relative to Ref. [2] , 
the angular distributions will be referenced to the median energy. The depth of 
the median energy is defined by 

Vo rVm 

aoiE{y))giy)dy^2 ao{E{y)) g{y) dy, (33) 

JO 

where ym is the median depth and the median energy is given by E„i = E{ym)- 
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Figure 9: The calculated correction factors for 48° (circles) 90° (squares), and 130° (diamonds) 
versus measurement number, for measurements 18-25. 



Note also that we have utilized the total ground state cross section in the def- 
inition. The deconvoluted differential cross sections, referenced to the median 
energies defined above, are calculated from the experimental target-averaged 



differential cross sections (defined by Eq. 29 ) via 

/dcr\ _ 1 /da\ 



\dn J f \dn/ ' ^^^^ 

where the correction factor is given by 

~ drr drr 

/=[A^.^(i^™)]-/ ^iE),iy)ay. (35) 

The theoretical form assumed for the differential cross section is the same as 
used in Subsec. [XSl 

'^(E) - 2^W,,iO). (36) 

The resulting correction factors are shown in Figs. [8] and [9] for three repre- 
sentative angles. For measurements 1-17 (below the peak of the narrow E2 
resonance), the correction factors are seen to differ from unity by at most 20% 
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and to depend very little upon the angle. For measurements 18-25, the cor- 
rection factors show much great departures from unity and significant angular 
dependence; these correction factors are also very sensitive to the assumed target 
depth profile. 

At this point, it is instructive to investigate how sensitive the deconvolution 
is to the assumed form of the cross section and depth profiles, both through the 
calculation of the median energy and the correction factor /. We will limit our 
consideration here to measurements 1-13, as measurements 14-17 are on the ris- 
ing slope of the narrow E2 resonance where taking the exact energy dependence 
of this resonance into account is clearly important, and it has already been noted 
that the measurements 18-25 are very sensitive to the assumed depth profile. 
To test the sensitivity to the assumed cross section, we have repeated the de- 
convolution assuming a constant S'-factor and cos 4>. In the case of the median 
energy, only measurements 1-2 show any sensitivity: the calculation with our 
model cross section yields median energies about 7-keV lower than the constant 
S'-factor calculation; note that this lowering corresponds to a a 4% increase in 
the calculated experimental S factors. The differences in the calculated cor- 
rection factors are at most 5%. We have also investigated the effects of a 20% 
increase in the scale of the depth profiles. This change leads to reductions in 
the median energy of up to 12 keV and up to 9% reductions in the correction 
factor. All these changes are within the tolerance implied by the uncertainties 
assigned to the effective energies given in column 2 of Table II of Ref . [5] . 

The deconvoluted experimental differential cross sections for measurements 
1-17 have been fitted using 



to determine <tei and aE2', cost/) has been fixed at the value given by the R- 
matrix parametrization shown in Fig. m The fit was determined by mini- 
mization, with the parameter errors indicating the range of parameter values 
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Figure 10: The El (upper panel) and E2 (lower panel) S factors. The values determined 
from the deconvoluted differential cross sections for measurements 1-17 are given as circles. 
The corresponding values given by Ref. [2j are given as squares; the error bars for these points 
are not shown but are essentially identical to the error bars on the corresponding circles. The 
i?-matrix parametrizations are given by the solid curves. 
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consistent with an increase in by an amount < 1 above the minimum when 
the remaining free parameters are varied. The uncertainties in N and the CM 
energy have been propagated after fitting. The resufts are shown as S'-factors in 
Figs. [I0| and [TT] and Table [T] For reference, the figures also show results given by 
Ref. from the analysis of the same raw data and the i?-matrix parametrization 
shown in Fig. |4j Our analysis differs from that given in Ref. |2j in three ways: 
the effects of CM motion have been taken into account, the median energy has 
been calculated using a more accurate representation of the energy dependence 
of the cross section, and the differential cross sections have been adjusted by the 
the correction factor given by Eq. |35j Another difference is that we arc fitting 
differential cross sections that have not been corrected for detector solid angle 
effects but rather include these corrections in our fitting function; this difference 
should not change the extracted S factors. The differences between our results 
and those given by Ref. [2] are generally quite small. In the case of the E2 S 
factor, the effects of CM motion and the deconvolution tend to cancel, leading to 
very small changes. Our results are also in generally good agreement with the 
i?-matrix parametrization, which indicates consistency between the data and 
our cross section model. Near the narrow E2 resonance, Ref. [2] provided two 
values of the energy (columns 2 and 3 of Table II of that reference); both values 
are plotted in Fig. [TT] We note that for measurements 14-17, both our results 
and those of Ref. [2] lie 3-10 keV to the left of the i?-matrix parametrization. 
This discrepancy is also seen in Fig. [6] The i?-niatrix parametrizaton should be 
well determined here by the known parameters of the resonance. One possibil- 
ity is that uncertainty in the beam energy contributes to this difference. In any 
case, the deconvolution procedure for measurements 14-17 must be viewed with 
some skepticism due to the lack of agreement between the resulting data points 
and our cross section model. 

In Fig. [12] we show the differential cross section obtained for measurement 7, 
for which we find median energy to be = 2.199 McV in the CM system. The 
two-parameter fit is shown as the solid curve; this fit utilized a fixed value of 
COS./) = 0.941 and yields asi = 24.3(5) nb and (7e2 = 0.8(1) nb with = 22.0. 
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E (MeV) 

cm. ^ ' 

Figure 11: The E\ (upper panel) and E2 (lower panel) S factors. The values determined from 
the deconvoluted differential cross sections for measurements 12-17 are given as circles. The 
corresponding values given by Ref. [2] using their energy method (I) are given as squares and 
for measurements 14-17 using their energy method (II) are as diamonds (shown only for the 
lower panel); the error bars for these points are not shown but are essentially identical to the 
error bars on the corresponding circles. The respective i?-matrix parametrizations are given 
by the solid curves. 
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4.0 




6 (degress) 

Figure 12: The deconvoluted angular distribution data for measurement 7. The sohd and 
dashed curves show the two-parameter and three-parameter fits that are discussed in more 
detail in the text. 

Assungao et al. [5] noted that three-parameter fits, where cos(f> was allowed 
to vary, tended to differ significantly from two-parameter fits for CM energies 
between 2.0 and 2.4 MeV. We find similar results. Such a three parameter 
fit for measurement 7 is also shown in Fig. [12] as the dashed curve; it yields 
asi = 21.9(8) nb, aE2 = 3.2(6) nb, and cos0 = 0.49(6) with = 8.3. The 
above uncertainties do not include the contribution from N which is common 
to all angles of a single angular distribution. It is seen that while both fits 
give essentially identical total ground-state cross sections (ctq = (^Ei + '^£2), 
the three-parameter fit gives a much larger aE2 ~ by a factor of 3.8 which 
corresponds to 3.6 standard deviations. We would like to emphasize that there 
is no known reason for the value for cos (f) assumed in the two-parameter fit to 
be incorrect (see discussion in Ref. [21 J. We have investigated the possible role 
of uncertainties in the deconvolution procedure on this finding. We conclude 
that it cannot be explained by this consideration, due to the lack of angle 
dependence of the convolution factor in this energy range (see Fig. Isl. The 
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most likely explanation is that angular distribution data contain unidentified 
systematic errors. Evidence for this explanation is provided by the ^ values 
from both the two-parameter and three-parameter fits to the data. Our 
values are very similar to those of Ref. |2j; as can be seen from Table I of 
Ref. [2], the values always exceed the number of degrees of freedom and 
often do so by a substantial amount. Another observation is that particular 
angles show systematic deviations from the fits. For example, the 130° data 
point lies above the fitted curve for 23 out of the 25 distributions in the case 
of the three-parameter fits (data near the narrow El resonance discussed in 
Subsec. [3^ were also considered). 

In Table [T] we also provide the total ground-state S factors 5*0, where the 
uncertainty is determined using two-parameter fits considering ctq and cteiIoex 
to be independent variables. The uncertainties in iV and the CM energy have 
again been propagated after fitting. In a simultaneous _R-matrix analysis of data 
from this experiment, we envision fitting fitting Sq and the angular distributions. 
The angular distributions would be considered to be un-normalized, with the 
normalization carried in the data, and an overall normalization uncertainty 
of 9% would also be assumed for the entire data set (inferred from Table 5.2 of 
Ref. |15j). This approach ensures that only independent data with independent 
uncertainties are fitted; this would not be the case if both the S'^i and Sei 
values were included in the fit. 

3. 7. Deconvolution of El data 

In Subsec. |3.5| it was noted that target-averaged El cross sections can be 
extracted independently of the E2 cross sections. This analysis will be pursued 
further here, as it is useful for the data near the narrow E2 resonance where 
the E2 contribution is extremely sensitive to the assumed target profile. In this 
case, we define the median energy using the El cross section: 



■V: 



<JEi{E{y))g{y) dy = 2 



(JEi{E{y)) g{y) dy, 



(38) 
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E 


Sei 


Se2 


So 




Sei 




(MeV) 


(keV-b) 


(keV-b) 


(keV-b) 


(MeV) 


(keV-b) 


1 


1.299 


20.8(85) 


4.3(54) 


25.1(73) 






2 


1.337 


14.5(59) 


11.0(51) 


25.5(69) 






3 


1.667 


20.5(34) 


4.8(12) 


25.3(39) 






4 


1.969 


26.5(39) 


3.6(7) 


30.1(43) 






5 


2.045 


30.1(44) 


3.9(7) 


34.1(49) 






6 


2.123 


41.5(60) 


2.5(5) 


44.0(63) 






7 


2.199 


57.6(82) 


2.0(4) 


59.5(84) 






8 


2.236 


61.1(86) 


2.2(4) 


63.3(89) 






9 


2.272 


79.0(110) 


2.9(5) 


81.9(114) 






10 


2.343 


71.0(99) 


1.5(3) 


72.6(101) 






11 


2.454 


43.8(31) 


2.3(4) 


46.1(32) 






12 


2.580 


19.7(24) 


2.0(7) 


21.7(25) 






13 


2.608 


16.8(15) 


1.8(4) 


18.6(16) 






14 


2.647 


10.9(12) 


7.1(10) 


18.0(16) 


2.645 


10.9(12) 


15 


2.655 


9.9(11) 


8.2(10) 


18.1(16) 


2.653 


10.0(11) 


16 


2.660 


10.1(15) 


11.4(15) 


21.5(20) 


2.656 


10.3(15) 


17 


2.666 


9.6(13) 


31.5(30) 


41.1(36) 


2.660 


10.2(14) 


18 










2.664 


8.3(28) 


19 










2.668 


10.8(29) 


20 










2.683 


7.7(19) 


21 










2.706 


7.4(13) 


22 










2.721 


4.8(9) 


23 










2.736 


4.8(9) 


24 










2.759 


3.8(5) 


25 










2.781 


4.0(5) 



Table 1: The deconvoluted 5-factor data. The first column supplies the measurement num- 
ber. The following four columns provide the El, E2, and total S-factor data from the two- 
parameter fits discussed in Subsec. |3.6| The final two columns provide the El 5-factor data 
from the three-parameter fits discussed in Subsec. |3.7| 



31 



where the median energy is given by Em = E(ym)- This procedure avoids 
introducing un-needed uncertainty through E2 contribution. The same target 



profiles are used as in Subsec. 3.6 The target-averaged El cross sections {(Tei) 



are obtained using three-parameter fits to Eq. (30). Dcconvolutcd El cross 
sections are obtained via 

(o-_Ei)cxp = -J— {(J El) (39) 
JEl 

where in this case the correction factor is given by 

fEi ^ [Ng<JEi{E„,)]-^ aEi{E)g{y)dy. (40) 
Jo 

The resulting El S factors for measurements 14-25 are given in Table [ij 

3.8. Electronic Files 

Three supplementary files accompany this paper on the publisher's website. 
These files contain the N^{0) data from Ref. [15], the Nj{6) data with the cor- 
rection for finite solid angle removed, and the N~^{0) data with the correction 
for finite solid angle removed and the deconvolution factor given by Eq. ( |35[ ) ap- 
plied, respectively. The latter data, when normalized by 0.624(47r)Q7V, become 
the deconvoluted experimental differential cross sections in nb/sr. 



4. Conclusions 

A general framework for deconvoluting the effects of energy averaging on 
charged-particle reaction measurements has been presented. There are many 
potentially correct approaches to the problem; the relative merits of some of 
them have been discussed. 

These deconvolution methods have been applied to recent ^^C(a, 7)^^0 mea- 
surements [5J[T5]. This analysis has clarified what can, and what cannot, be 
explained by energy convolution effects. Very significant effects are found at 
and above the narrow 2.68-MeV resonance. Below this resonance, we find that 
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considerations of energy deconvolution, as well as special relativity, lead to rel- 
atively small changes in the extracted S factors. We have also extracted decon- 
voluted differential cross sections from this experiment. We expect that these 
results will be useful for future analyses, including a simultaneous i?-matrix of 
the angular distributions. 
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